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Abstract. Continuing the lines developed in [2T], in this paper we study nonhomogeneous wavelet 
systems in high dimensions. It is of interest to study a wavelet system with a minimum number 
of generators. It has been showed by X. Dai, D. R. Larson, and D. M. Speegle in JTTJ that 
for any d x d real- valued expansive matrix M, a homogeneous orthonormal M -wavelet basis can 
be generated by a single wavelet function. On the other hand, it has been demonstrated in 
[2"T] that nonhomogeneous wavelet systems, though much less studied in the literature, play a 
fundamental role in wavelet analysis and naturally link many aspects of wavelet analysis together. 
In this paper, we arc interested in nonhomogeneous wavelet systems in high dimensions with a 
minimum number of generators. As we shall see in this paper, a nonhomogeneous wavelet system 
naturally leads to a homogeneous wavelet system with almost all properties preserved. We also 
show that a nonredundant nonhomogeneous wavelet system is naturally connected to refinable 
structures and has a fixed number of wavelet generators. Consequently, it is often impossible for 
a nonhomogeneous orthonormal wavelet basis to have a single wavelet generator. However, for 
redundant nonhomogeneous wavelet systems, we show that for any d x d real- valued expansive 
matrix M, we can always construct a nonhomogeneous smooth tight M-wavelet frame in L 2 (R d ) 
with a single wavelet generator whose Fourier transform is a compactly supported C°° function. 
Moreover, such nonhomogeneous tight wavelet frames are associated with filter banks and can be 
modified to achieve directionality in high dimensions. Our analysis of nonhomogeneous wavelet 
systems employs a notion of frequency-based nonhomogeneous wavelet systems in the distribution 
space. Such a notion allows us to completely separate the perfect reconstruction property of a 
wavelet system from its stability in various function spaces. 



For a function / : M. d — > C and a d x d real-valued invertible matrix U, throughout the paper 
we shall adopt the following notation: 



fu;wAx):=\detU\ 1/2 e- m - Ux f(Ux-k) and := / w , x, k, n G R d , (1.1) 



where i denotes the imaginary unit. Now we recall the definition of a homogeneous wavelet system, 
which is closely related to the discretization of a continuous wavelet transform. Let M be a d x d 
real- valued invertible matrix and let ^ be a subset of square integrable functions in L2(M. d ). The 
following homogeneous M-wavelet system 



has been extensively studied in the function space L2(R d ) in wavelet analysis, often with M being 
an integer expansive matrix. Here we say that M is an expansive matrix if all its eigenvalues have 
modulus greater than one. The elements in \l/ of fll.2j) are called wavelet functions or wavelet 
generators. It is important to point out here that the elements in a set S of generators in this 
paper are not necessarily distinct and S may be an infinite set. The notation h G S in a summation 
means that h visits every element (with multiplicity) in S once and only once. For a set S, we 

2000 Mathematics Subject Classification. 42C40, 42C15, 42B05. 

Key words and phrases. Wavelet systems with a minimum number of generators, nonhomogeneous wavelet 
systems, homogeneous wavelet systems, frequency-based nonhomogeneous dual wavelet frames, distribution space, 
Fourier transform, nonstationary wavelet systems, real- valued dilation matrices. 

Research supported in part by NSERC Canada under Grant RGP 228051. February 11, 2010. 



1. Introduction and Motivations 



WS(vI/j: {c Ml:< : J G Z, k G 




l 



2 



BIN HAN 



shall use #5 to denote its cardinality counting multiplicity. For example, for \l/ = {ijj 1 , . . . ,if> s }, 
its cardinality is s and s could be oo, all the functions if) 1 , . . . , if) s are not necessarily distinct, 
and if) G ^ in (11 .2p simply means ^ = . . . For some references on homogeneous wavelet 
systems, see [U|2], [5] [30]. 

It is of interest to study a homogeneous wavelet system with a minimum number of generators. 
It has been showed in the interesting paper of Dai, Larson, and Speegle [HI Corollary 1] that 
for any real- valued expansive matrix M, there exists a (Lebesgue) measurable subset E of IR d 
such that the homogeneous M-wavelet system WSdV'}) is an orthonormal basis in L2(M d ), where 
if, = xe, the characteristic function of E. Such a measurable set E is called a dilation-M wavelet 
set in [10J [11] . For more discussion on wavelet sets and their constructions, see [2j [TO], [TTJ [25] and 
references therein. Since an orthonormal wavelet basis is a special case of a tight (or Parseval) 
wavelet frame, such wavelet sets are automatically framelet sets which yield homogeneous tight 
wavelet frames in L 2 (M. d ) with a single wavelet generator. For M = 21 ^ where Id denotes the 
d x d identity matrix, it has been showed in [ISJ Theorem 3.8] that one can always construct a 
homogeneous smooth tight M-wavelet frame WS({?/>}) with a single generator such that if) is a 
compactly supported function in C°°(lR a! ). The Fourier transform used in this paper for / G Li(IR d ) 
is defined to be /(£) = J Rd f{x)e~ lx '^dx, £ G R d and can be naturally extended to square integrable 
functions and tempered distributions. Through the study of fast wavelet transforms and their 
underlying wavelet systems, it has been recently noticed in [2TJ that a nonhomogeneous wavelet 
system, though much less studied in the literature, plays a fundamental role in understanding many 
aspects of wavelet analysis and its applications. In this paper, we shall continue the lines developed 
in [21] for dimension one to investigate nonhomogeneous wavelet systems in high dimensions. Let 
$ and \1/ be subsets of L2(M. d ). A nonhomogeneous M-wavelet system is defined to be 

WSj(^ *) := {0 MJ;k : k G Z d , G <£} U {^k : J > I k G Z d , if) G (1.3) 

where J is an integer representing the coarsest scale level. In this paper, we are particularly 
interested in nonhomogeneous wavelet systems with a minimum number of generators, that is, 
the smallest possible cardinalities and of generators in (jl.3p . In contrast to homogeneous 
wavelet systems, we shall see in this paper that there is an intrinsic difference between nonredun- 
dant and redundant nonhomogeneous wavelet systems. To have some rough ideas about our results 
in this paper, let us present here two results on nonredundant and redundant nonhomogeneous 
wavelet systems with a minimum number of generators. 

For nonhomogeneous orthonormal wavelet bases, we have the following result, which is a special 
case of Theorem [71 

Theorem 1. Let M be adxd integer invertible matrix. Let $ = {(f) 1 , . . . , cff} and = {if) 1 , . . . , ip s } 
be subsets o/L 2 (M d ) with r,s G NU {0}. Suppose that the nonhomogeneous M-wavelet system 
WSj($; ^) is an orthonormal basis o/L 2 (M d ) for some integer J. Then s = r(| det M| — 1) and 
there exist an r x r matrix a and an s x r matrix b of 2n -periodic measurable functions in L 2 (T d ) 
such that 

0(M T O = a(O0(O and ip(M T ^) = b(£)0(£), a.e.£eR d , (1.4) 

where (f> := [(f) 1 , ... , (p r } T and if) :— [ip 1 , . . . , if) s ] T . If in addition M is an expansive matrix, then the 
homogeneous M-wavelet system WS(\I / ) is an orthonormal basis of L2(M. d ). 

By Theorem [TJ for a nonhomogeneous orthonormal M-wavelet basis in Zv2(M d ), the smallest 
possible numbers of generators are = 1 and = 1, for which the integer matrix M must 
satisfy | det M| = 2. 

In contrast to nonredundant nonhomogeneous wavelet systems in Theorem [TJ, for redundant 
nonhomogeneous wavelet systems such as tight wavelet frames (and consequently dual wavelet 
frames), we have the following result, whose proof and construction are given in Section 4. 
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Theorem 2. Let M be a d x d real-valued expansive matrix. Then there exist two real-valued 
functions <j), ip in the Schwarz class such that 

(i) WSj({0}; {ip}) is a nonhomogeneous tight M -wavelet frame in 1^(18^) for all integers J: 

oo 

\\f\\l m = E \(f>M\ 3 + E E M'Ml* V / G L 2 (R d ); (1.5) 
kez d j=J kez d 

(ii) (j) and ijj are compactly supported C°° even functions; 

(iii) if) vanishes in a neighborhood of the origin, i.e., ip has arbitrarily high vanishing moments; 

(iv) there exist 27rZ d '-periodic measurable functions a^, b^, k e 7L d in C°°(T d ) such that 

e- lkMT ^(M T O=a k (O0(O and e -* kMT ^(M T = b k (OHO, HeR d ,keZ d . (1.6) 
Moreover, WS({ip}) is a homogeneous tight M -wavelet frame in L 2 (R d ) satisfying 

1 -EEl^'^W)! 2 V/GL 2 (M d ). (1.7) 



L 2 

The property in (11.61) is called the refinable structure of and ^ in this paper. We shall see in 
Section 4 that the nonhomogeneous tight M-wavelet frames in Theorem [2] can be easily modified 
to achieve directionality in high dimensions by using nonstationary tight wavelet frames. 

The structure of the paper is as follows. In Section 2, we shall explore the connections between 
nonhomogeneous and homogeneous wavelet systems in L 2 (IR d ). We shall see in Section 2 that any 
given nonhomogeneous wavelet system will yield a homogeneous wavelet system and a sequence 
of nonhomogeneous wavelet systems with almost all properties preserved. For a nonredundant 
nonhomogeneous wavelet system, we shall see in Section 2 that it has a natural connection to 
refinable structures. In Section 3, we shall introduce and characterize a pair of frequency-based 
nonhomogeneous (and more generally, nonstationary) dual wavelet frames in the distribution 
space. We shall see in Section 3 that this notion allows us to completely separate the perfect 
reconstruction property of a wavelet system from its stability in various function spaces. As an 
application, for any d x d real- valued invertible matrix M, we obtain a complete characterization 
of a pair of nonhomogeneous (as well as nonstationary) dual wavelet frames in /^(M^). This 
also yields a characterization of nonstationary tight wavelet frames in L2(M d ). Based on the 
results in Section 3, we shall prove Theorem [2] in Section 4. Moreover, we shall show that such 
nonhomogeneous tight wavelet frames in Theorem [2] are associated with filter banks and can be 
easily modified to achieve directionality in high dimensions. 



2. Nonhomogeneous and Homogeneous Wavelet Systems in L 2 ( 



i)d\ 



In this section, we shall study the connections of a nonhomogeneous wavelet system in L 2 (IR d ) 
to a homogeneous wavelet system in L 2 (M. d ). To do so, we need the following auxiliary result. 

Lemma 3. Let M be a dxd real-valued expansive matrix and $ be a (not necessarily finite) subset 
of L 2 (R d ) such that ^^ e<J) 1 1 1 1 x, 3 < 00 ■ Suppose that there exists a positive constant C such 
that E^Ekez' !(/,</»(• - k))| 2 ^ C\\f\\l md) for all f E L 2 (R d ). Then 

hm VV|(/,^ k )| 2 =fl V/GL 2 (M d ). (2.1) 

Proof. We use a similar argument as in [2TI Lemma 3]. We first prove that (12.11) holds with f = Xe, 
the characteristic function of a bounded measurable set E. By calculation, for / = xe, we have 

\(f,(f) W . k )\ 2 = |detMp f <p(Wx-k)dx 2 ^ | det M|- J ( I \(f)(x)\dx) 2 <: \E\ [ \4>(x)\ 2 dx, 

JE ^JMiE-k ' JWE-k 
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where \E\ denotes the Lebesgue measure of E. Since E is bounded and M is an expansive matrix, 
we see that M^E — k, k G Z d are mutually disjoint as j — > —oo. Therefore, we have 

EE^^^^EE / \<Kx)\*dx = \E\ [ h(x)dx, 

#E$ keZ d keZd JWE-k Ju kezd (WE-k) 

where h(x) := 5^ e $ |0(x)| 2 . By our assumption on $, we have h G Li(R d ). Since ft, G Li(lR d ), 
for arbitrary e > 0, there exist positive constants N and c such that J{ yl zM d -\y\>N} h(x)dx < e and 
j K h(x)dx < e for every measurable set K with \K\ < c. Since M is an expansive matrix, there 
exists an integer J such that for all j ^ J, the Lebesgue measure of the set K N := {y G M, d : \y\ ^ 
iV} fl (U keZ[ ^(M : '_E , — k)) is less than c. Consequently, we deduce that 

^2 E l(/'^Mi;k)| 2 ^ \E\ / h(x)dx^\E\ h(x)dx+\E\ h(x)dx^2\E\e 

for all j ^ J. Therefore, we see that ( 12.1 p holds for / = xe- Consequently, ( 12. ip holds for 
any function / which is a finite linear combination of characteristic functions of some bounded 
measurable sets. 

Define operators Pj : L 2 (R d ) -> / 2 (Z x $) by Pjf := {(/, 0Mi;k)}k£Z«*,06#- Since (/>0iW;k) = 
(/ M - J;O ,0(- — k)), we can easily deduce that 

ii^/ii 2 : = E E k/><w>i 2 = E E i</m-w(- - k)>r < cii/ M -, ;0 in 2(Rd) = cu/u^, 

for all j G Z. For any / G L2(M d ) and arbitrary e > 0, there exists a function g, which is a finite 
linear combination of characteristic functions of bounded measurable sets, such that ||/— <?||wRd) ^ 
e. Since we proved that lim^.oo ||Pj^|| = 0, there exists J G Z such that \\Pjg\\ ^ £ for all j ^ J. 
Consequently, we have 

11^711 < ll^ll + II W -g)\\^£ + y/c\\f- g\\ L2(Rd) ^e(i + Vc), v j ^ J. 

Therefore, we conclude that lim 7 _ i ,_ 00 ||Pj/|| = 0. That is, (12. ip holds. □ 

We mention that ( 12. ip in Lemma |3] holds for a more general d x d real- valued (not necessarily 
expansive) matrix. For example, by a slightly modified proof of Lemma |3l (12. ip holds if M is a 
d x d real-valued matrix such that 

| det M | > 1 and limsup sup ||M J x|| < oo. (2.2) 

j— >-oo ||a;||$l 

Without requiring that M be an expansive matrix, all the results in this section still hold provided 
that (12.11) holds even with lim in ( 12. ip being replaced by liminf. 

Proposition 4. Let M be a dx d real-valued invertible matrix. Let $ and \1/ be subsets of L2(M d ). 
Suppose that WS j($; ty) is a nonhomogeneous M-wavelet frame in /^(IR^) for some integer J, that 
is, there exist positive constants C\ and G% such that 

oo 

Ci\\f\\i 2{Rd) < E E k/.wi 2 + EEE i(/^M, ; k)i 2 < c 2 \\f\\i 2(Rd) , 



V/6l 2 (l d ). 

Then (I2.3P holds for all integers J, in other words, WSj(<3>; is a nonhomogeneous M-wavelet 
frame with the same lower and upper frame bounds for all integers J. If in addition M is an 
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expansive matrix and X^ e $ ll^llz^pf^) < 00 > then the homogeneous M-wavelet system WS(\1/) is 
also a frame in L 2 (M. d ) with the same lower and upper frame bounds: 

Ci\\f\\l 2m ^ EE E K/>^ ;k >r < C 2 \\f\\l m , V/ G L 2 (R d ). (2.4) 

jez vg* kez d 

Proof. By the following simple fact that for f,g G L 2 (IR d ) and invertible <i x d matrices U and V^, 

(fu;k, n , gu-,k, n } = (/, #) and (/y ;0 ,o, ^a ; k,n) = (/, ^c/y-i ; k,n), k, n G Z d , (2.5) 

it is straightforward to see that if (12 .3p holds for one integer J, then (I2.3P holds for all integers J. 
Next we prove that (12 .4p holds. Since (I2.3P holds for all integers J, it is easy to see that 

00 

EEEK/'Wr^^n/iiL^) 

j=j if)£-t> kez d 

and 

Ci\\f\\i 2i ^ -EE ia^ ;k )i 2 < E E E K/'Wi 2 

for all integers J G Z and / G L 2 (IR d ). Taking J — > — 00 in the above two inequalities and using 
( 12. ip for the second inequality, we deduce that (12.4j) holds. □ 

The best possible constants C\ and C 2 in ( 12. 3 p are called the lower frame bound and the upper 
frame bound of WSj($; \I/), respectively. 

Let M be a d x d real-valued invertible matrix. Let 

$ = {0\...,<f}, * = {V>V..,^} and $ = {0 1 ,...,0 r }, * = tfV..,^} (2.6) 

be subsets of L 2 (R d ), where r,sG NU{0, +00}. Let WSj($; tf) be defined in flTSJ) and WSj(<§; ^) 
be defined similarly. We say that the pair (WSj($; \I/), WSj($; \&)) is a pair of nonhomogeneous 
dual M-wavelet frames in L 2 (R d ) if each of WSj($; \1/) and WSj($; \l/) is a nonhomogeneous M- 
wavelet frame in L 2 (R d ) and the following identity holds 

r 00 s 

(f>9) = E E(/'A^ft J M) + EEE</'ikHft ik -ff>> V/,0 G L 2 (R d ). (2.7) 

For pairs of nonhomogeneous dual wavelet frames in L 2 (M d ), we have the following result. 

Proposition 5. Let M be a d x d real-valued invertible matrix. Let $, $, \l/ zn ( 12. 6 p &e subsets 
of L 2 (R d ). Suppose that (WSj(<3>; \l/), WSj($; zs a pazr of nonhomogeneous dual M-wavelet 
frames in L 2 (R d ) for some integer J. Then it is a pair of nonhomogeneous dual M-wavelet frames 
in L 2 (R d ) for all integers J. If in addition M is an expansive matrix and X^e$ II ^11 z, 2 (K d ) < 00 ' 
then the pair (WS(\I/), WS(\I/)) is a pair of homogeneous dual M-wavelet frames in L 2 (R d ), that 
is, each of WS(\I / ) and WS(\I/) is a frame in L 2 (R d ) and the following identity holds: 

s 

if' 9) = E E E (/' V4i,k> *;k> 9), V f,g G L 2 (M d ). (2.8) 

Proo/. By (|2~3|l . we see that (J22D holds for all integers J. By PropositionH (WSj($; WSj($; *)) 
is a pair of nonhomogeneous dual M-wavelet frames in L 2 (M d ) for all integers J. For a real- valued 
expansive matrix M, by Proposition H] again, each of WS(\I/) and WS(\I/) is a frame in L 2 (R d ). By 
Lemma EJ the identity in ( 12. 8ft can be proved similarly as in the proof of Proposition HI □ 



6 



BIN HAN 



Next, we discuss nonredundant nonhomogeneous wavelet systems. Let $ and \& be subsets of 
L 2 (M. d ). We say that WSj($; \J r ) is a nonhomogeneous Riesz M-wavelet basis in L 2 (M. d ) if the linear 
span of WSj($; \E r ) is dense in L 2 (M. d ) and there exist two positive constants C3 and C4 such that 



Ci { Y Y \ v K<p\ 2 + Yj Y Y \ w iM\j <|| Y Y v K^MJ;k + Y Y Y w J;K^Mi;k 

00 

^ Ci {Y Y M 2 + YY Y K;M-I 2 



Li 



(2.J 



for all finitely supported sequences {i>k,</>}kez d ,4>e$ an d {wjik^jj^j.kez^Vie*) where the best possible 
constants C3 and C4 are called the lower Riesz bound and the upper Riesz bound of WSj($; \&), 
respectively. It is well known that WSj(<3>; is a Riesz basis of L 2 (M d ) if and only if it is a frame 
of L 2 (M. d ) and it is l 2 -linearly independent, that is, if 

00 

Y Y v ^MJ;k + YYY w i;M^ ; k = (2.10) 
06* kez d j=J V'e* kez d 



in L 2 (M ) for square summable sequences {fk,</>}kez d ,</>e<i> an d {^k^jj^kez^e*' then = and 
Wj . M = for all j ^ J, k G Z d and G $, ip G In fact, define W : Z 2 (Z d x$)xl 2 (ZxZ d xf) 4 
L 2 (IR d ) such that W maps ({fk,0}kez d ,</>e*; { w j;k,^}j3 ! ,/,kez d ,v>e'i') to the function on the left-hand 
side of ( l2~T0i) . Then WSj($; \I>) is a Riesz basis in L 2 (R d ) if and only if W is a well-defined 
bounded and invertible operator. Also, WSj($; is a frame in L 2 (M d ) if and only if W is a 
well-defined bounded onto operator with the range of its adjoint operator W* being closed. Now 
it is straightforward to see that WSj($; is a Riesz basis of L 2 {R d ) if and only if it is a frame 
of L 2 (M. d ) and it is / 2 -linearly independent. Thus, a nonhomogeneous Riesz M-wavelet basis is 
nonredundant. Note that a nonhomogeneous orthonormal M-wavelet basis is a particular case of 
a nonhomogeneous Riesz M-wavelet basis. 

For nonhomogeneous Riesz wavelet bases in L 2 (M, d ), we have the following result. 

Theorem 6. Let M be a d x d real-valued invertible matrix. Let $ and ^ be subsets of L 2 (M. d ). 
Suppose that WSj($; is a nonhomogeneous Riesz M-wavelet basis in L 2 (M. d ) satisfying (12. 9p for 
some integer J. Then (12. 9p holds for all integers J and WSj($; is a nonhomogeneous Riesz 
M-wavelet basis in L 2 (M. d ) with the same lower and upper Riesz bounds for all integers J. If M 
is a d x d real-valued expansive matrix and X^e* ll^llL 2 (iR d ) < °°> then WS(\1/) is a homogeneous 
Riesz M-wavelet basis in L 2 (M. d ) with the same lower and upper Riesz bounds, that is, the linear 
span o/WS(\l/) is dense in L 2 (JBL d ) and 

11 2 

C3 YYY k^i 2 < || Y Y Y „ < C4 Y Y Y k^i 2 ( 2 - n ) 

j& ve* kez d jez i/>e* kez d 2 jez kez d 

for all finitely supported sequences {w,j;k,v>}jez,kez <i W>e*- 

Proo/. Since (El]) holds, by (V>Mi ; k)M» ; o = ^Mi+^k and ||/iW;o|li 2 ( K d) = ||/Hwr<i)» it is eas y to veri fy 



that (E3D holds for all integers J. Note that WS J+n ($; *) = {/ M » ; o : / e WSj($;^)} for all 
integers J and n. It is also easy to deduce that if the linear span of WSj($; ^) is dense in L 2 (M. d ) 
for one integer J, then it is also dense in L 2 (M. d ) for all integers J. Consequently, WSj(<I>; is 
a nonhomogeneous Riesz M-wavelet basis in L 2 (R d ) for all integers J with the same lower and 
upper Riesz bounds. 

Since (12. 9p holds for all integers J, setting = in (I2.9p . we can easily deduce that (12. lip 
holds for all finitely supported sequences {^j ; k,i/j}j£Z,kez d ,</>e*- Since WSj(<I>; is a Riesz basis in 
L 2 (M. d ), we see that WSj($; ^) must be a frame in L 2 (M. d ). Since M is an expansive matrix and 
X^e* 1 1 ^ 1 1 2,2 (M d ) < 00 ' ^ Proposition HI WS(^) is a frame in L 2 (M d ). Therefore, by (12. 4p . we see 
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that / _L WS(\I/) implies / = 0. Thus, we conclude that the linear span of WS(\I/) is dense in 
L 2 (R d ). Therefore, WS(tf) is a Riesz basis of L 2 (R d ). □ 

Let 5 denote the Dirac sequence such that 5(0) = 1 and S(k) = for all k ^ 0. Let <E>, $, \I/ 
in fLT6|) be subsets of L 2 (IR d ). We say that the pair (WSj($; WSj(l>; *)) is a pazr o/ non- 
homogeneous biorthogonal M-wavelet bases in L 2 (M d ) if each of WSj($; and WSj($; "3*) is a 
nonhomogeneous Riesz M-wavelet basis in L 2 (M. d ) and the following biorthogonality relation holds: 

(A^*M^>=*(k-W-0, » ;k ,& ;k ') = 0, (0^<y ;k ^) = O, (2.12) 

» ;k ,</ ;k ,> = 5(k-k')5(n-n')5(j-/), (2.13) 

for all k, k' G Z d , j, j' G Z such that j, f ^ J, £,£' = 1, . . . ,r, and n, n' = 1, . . . , s. 

It is a standard result, which can be easily proved by the same argument after (I2.10p . that a pair 
(WSj(<3>; \1/), WSj($; \I/)) is a pair of nonhomogeneous biorthogonal M-wavelet bases in L 2 (R d ) if 
and only if it is a pair of nonhomogeneous dual M-wavelet frames in L 2 (M d ) and the biorthogonality 
conditions in (I2.12p and (12.131) are satisfied. For pairs of nonhomogeneous biorthogonal M-wavelet 
bases in L 2 (M. d ), we have the following result which includes Theorem [1] as a special case: 

Theorem 7. Let M be a d x d real-valued invertible matrix. Let $, $, in (12. 6p be finite 
subsets of L 2 (M> d ). Suppose that the pair (WSj($; \1/), WSj($; \1/)) zs a pair of nonhomogeneous 
biorthogonal M-wavelet bases in L 2 (M. d ) for some integer J. Then it is a pair of nonhomogeneous 
biorthogonal M-wavelet bases in L 2 {U. d ) for all integers J. Moreover, there exist r x r matrices 
a^, a^ and s x r matrices b^, b^, k G Z d of 2^7^ -periodic measurable functions in L 2 (T d ) such that 

e-* k - MT «0(M T O = a k (O0(O and e" lk MT ^(M r O = b k (£)0(£), a.e. £ G M d , k G Z d , (2.14) 
e- k ' MT «0(M T O = a k (£)|(£) and e" tk MT ^(M T = b k (£)0(£), a.e. £ G M d , k G Z d , (2.15) 

0:=[0 1 ,...,0r, V> := [V>\ • • • ,^F, 0:=[0 1 ,-..,0T, ^ := 1 , • • • , ■ (2-16) 
// M is a d x d integer invertible matrix with cIm := | det M| ; then s = r(d|\/i — 1) and 



(2.18) 



P[ao,bo](O^P[ao,b ](0 = /rd M> a - e - £ ^ (2-17) 

where I r ^ M denotes the (rdu) X (rdu) identity matrix and 

P (f=\ = T a o(£ + 2tto;o) a (£ + 27rwi) ■•• a (£ + 27ra; dM _i) 
nao,bo]W- [b (£ + 27ra;o) ^(£ + 2^) •■■ b (£ + 27ra; dM _i; 

and {c^o, . . . , a;d M -i} := [(M T ) _1 Z d ] fl [0, l) d . // M is a d x d real-valued expansive matrix, then 
(WS(\I/), WS( 1 P)) a pazr of homogeneous biorthogonal M-wavelet bases in L 2 (R d ), that is, each 
of WS(\I/) and WS(\I/) a homogeneous Riesz M-wavelet basis in L 2 (M. d ) and ( I2.13P holds for all 
k, k' G Z d , j, f G Z ; and n, n' = 1, . . . , s. 

Proof. By Theorem [6j each of WSj($; and WSj($; ty) is a nonhomogeneous Riesz M-wavelet 
basis in L 2 (R d ) for all integers J. By (12. 5p . we see that (I2.12p and ( I2.13P must hold for all integers 
J. Thus, (WSj($; \1/), WSj($; \l/)) is a pair of nonhomogeneous biorthogonal M-wavelet bases in 
L 2 (R d ) for all integers J. 

To prove ( I2.14p and ( I2.15p . let us consider the representations of 0^_i. k , ^^_i. k , k G Z d , £ = 
1, . . . , r and n = 1, . . . , s under the Riesz basis WS ($; Noting that ( 12TT2"]) and'()27l3]) hold for 
all integers J, we deduce that 

r r 

^M-i;k = X X (^M-l; k) 4>t, m )<f>t, m and VB|-i;k = X) 2 Wl-^lo ^l;m)<Cm ( 2 -19) 
<'=1 mGZ d «'=1 m€Z d 
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with all the coefficient sequences being square summable and the series in (I2.19p converging in 
L 2 (R d ). For k £ Z d , define 



(2.20) 

MOW' := | det Ml" 1 / 2 £ (fr ik ,t)c- § , 

where [a^]^/ denotes the (£, £')-entry of the matrix a^. It is evident that both a k and b|< are 
matrices of 27rZ d -periodic measurable functions in L 2 (T d ). Taking Fourier transform on both 

sides of the equations in ( 12 . 19[) and noting 0^-1. k (0 = I det M| 1//2 e _:(k ' MT ^0(lv1 T ^), we conclude 

that (I2.14p holds. Using the same argument by switching the roles of <f), ip with <f>, ijj, we see that 
( 12T51) holds with 

[MOW := | det Mr 1 / 2 X)#M-;k.<;m>e- im - «, 

m&d (2.21) 
MOW' := I det Ml" 1 / 2 (A-i ; k.<;m>e- <m - 5 - 

mGZ d 

For an integer invertible matrix M, we denote by {70, . . . , 7d M -i} := {Mx : x £ (M _1 Z d ) fl 
[0, l) d }. That is, {70, . . . , 7d M _i} is a complete set of representatives of the distinct cosets in the 
quotient group Z d /(MZ d ). For a (0 = Ekez" a (k)e~ ik <, we have a (0 = En=o* a , 7n (M T 0e-^ 
with a 0)7n (0 := Ekez d °o(Mk + 7 n )e~* k ''». Now it is easy to check that 

[a (e + 27ru;o),...,ao(e + 27ra; dM _ 1 )] = [a , 70 (M T £), . . . , ao^.^M^MO, (2.22) 

where 

U(0 := (e- i(5+27ra)m) ^/,)o< w> ^dM-i and U(0U(0 T = d M W 

Denote 



/i := [0\...,0 r ,V\---,^ s ] T and 77 := [0^ ;7o , . . . , <^, 7o , . . . , 4>u. yiu _ v ■ ■ • , 0M, 7 d M -i lT 



Noting that 0m ;7 „(O = d M 1/2 e- i7 " N€ <//(N£) with N := (M T )W we deduce from (l2~T9l) and (1^221) 
that 

MO = d^P^^KNOU-^NO^O- (2-23) 
Representing the entries in the vector rj under the Riesz basis WSo($; ^), we have 

<4 7 „ = EE ^M; 7 „, 0l;-kX;-k + EE (^; 7n) ^l-kX^k. (2-24) 

£'=1 kez d n'=i kez d 

Using OMU) and <^M>, by (<^ ;7 „, <^._ k ) = (^ M;Mk+7n , ¥), we deduce that 

m = dM /2 U r W T P [iio ,bo](N0 T M0 = dM 1/2 U(N0P [io ,6 o] (N0 r ^(0- (2-25) 

Note that both — k), k £ Z d and rj(- — k), k £ Z d are two Riesz bases for the same subspace. 
Now it follows from (12351) and (T2T25]) that 

-T 



P[ao,bo](N£)P[ao,bo]( N = L r + s and (N£) P [a0)bo] (N£) = Jhm • 

The above identities hold if and only if r + s = rdf,/\ (that is, s = r (cJm — 1)) and (12.171) holds. 

For a real- valued expansive matrix M, by Theorem [61 each of WS(\I/) and WS(\I/) is a Riesz 
basis in L 2 (R ci ). Since ( 12. 13ft holds for all integers j and f, we conclude that (WS(^), WS(\E f )) is 
a pair of homogeneous biorthogonal M-wavelet bases in L2(M d ). □ 
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Using an argument in [19) Lemma 1], we explore in the following result the connections between 
nonhomogeneous Riesz wavelet bases and refinable structure. 

Theorem 8. Let M be a d x d integer invertible matrix. Let $ and \1/ in (12. 6p be finite subsets of 
L 2 (IR d ). Suppose that WSj(<&; \P) is a nonhomogeneous Riesz M-wavelet basis in L 2 (IR d ) for some 
integer J. Then the following statements are equivalent to each other: 

(i) there exist subsets $ and ^ in ( T276T) of L 2 (M d ) snca that the pair (WSj($; WSj(4>; §)) 
zs a pazr of nonhomogeneous biorthogonal M-wavelet bases in L 2 (IR d ); 

(ii) there exist an r x r matrix a and an s x r matrix b o/ 27rZ d -periodic measurable functions 
in L 2 (T d ) such that (11. 4p ao/ds wift and ^ feeing defined in (12. 16)) . 

Proof. (i)=^(ii) is a direct consequence of Theorem)?) We now prove (ii) (i) . Since WSj($; is 
a Riesz basis in L 2 (R d ), it has a unique dual Riesz basis {0 £ ' ;J,k ' : k' G Z d , = 1, . . . , r}U{^ n ' ;: ''' k/ : 
k' G Z d , j' ^ J,n' — 1, . . . , s} such that (12.121) and (12. 13)) hold with 0^j. k , and V>J}|j'-k' being replaced 
by <p e '' J ' k ' and t/j n '' : >'' k \ respectively. Define 

^ : =#m-3o. ^ = l,---,r and ^' := ^'i J j° , n' = 1, . . . , s. (2.26) 
Now we prove that (12. 12j) and (12. 13j) must hold. In fact, by the definition in (12.26)) . we have 

(A^k.#M^> = (0M-k-k'>0M-o) = (0W-k^ V '°> = 5(£-l')5(k-W) 
and for j ^ J, noting that IVP _J Z d C Z d , we have 

(^Mi;k; ^M J ;k') = (^W ;k-IVp- J k' > ^M J ; o) = (^M* jk-lW" J k' > ^ ' J ' ) = 0- 

We also observe that 

(0M';k.$iV> = (^M-k-k'^MVo) = (0M-k-k^ nV '°) = 0. (2.27) 

We now prove the rest of (I2.12p and (12.131) by a similar argument as in [191 Lemma 1]. For 
f > J, by (11. 4p . we see that 0^j. k is an / 2 -linear combination of 0^/. , m G Z d and L = 1, . . . , r. 

By (|22ID, we have (^. m ,< y . k ,> = (^ J;m ,^V> = °- Therefore,' (/ M . J;k , ^,. k ,) = 0. Hence, 
flZTD is verified. 

To prove (12. 13j) . we consider two cases. If j ^ j' ^ J, then 

= (^ + , ;k _ M y-i' k '^ nV,0 > = *0'-mk-lO- 

If j' > j ^ J, then by (11 ,4p . ^jvp'-k ^ s an ^-linear combination of </>J^j'- , m G Z d and L = 1, . . . , r. 

By (J22ZD, we have (^. m ,^. k ,) = (^ ;m A ; k'> = 0. Therefore,' (V^ ;k ,VftV. k '> = °- Hence > 
(12TT31) is verified. 

In conclusion, by the uniqueness of a dual Riesz basis of WSj($; \&), we proved that (fi e M j. k , = 
0£'-,J,k' anc j ^"'_._ k/ — Tpn';j,k> f Qr a jj j/ ^ ^ k' g a nd f = 1, . . . , r, n' = 1, . . . , s. Therefore, 

the pair (WSj($; \l/), WSj($; \&)) is a pair of nonhomogeneous biorthogonal M-wavelet bases in 
L 2 (IR d ). This completes the proof of (ii)=>(i). □ 

We conclude this section by some remarks. Firstly, by the results in this section, we see that a 
nonhomogeneous wavelet system WSj($; in L 2 (IR d ) for a given integer J (that is, at the coarsest 
scale level J) naturally leads to a sequence of nonhomogeneous wavelet systems WSj($;^/) for 
all integers j G Z while preserving almost all the properties of the system at j = J. This is 
a fundamental property in wavelet analysis. In fact, a one- level fast wavelet transform is just 
a transform between the two sets of wavelet coefficients of a given function represented under 
two nonhomogeneous wavelet systems at two consecutive scale levels. Naturally, for a multi-level 
wavelet transform, there is an underlying sequence of wavelet systems at every scale level, instead 
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of just one single wavelet system. For a given nonhomogeneous (stationary) wavelet system, since 
it naturally produces a sequence of nonhomogeneous wavelet systems, this allows us to study only 
one nonhomogeneous wavelet system instead of a sequence of nonhomogeneous wavelet systems. 
This desirable property of nonhomogeneous wavelet systems is not shared by homogeneous wavelet 
systems. Secondly, by the results in this section, for any d x d real-valued expansive matrix 
M, we see that a nonhomogeneous wavelet system WSj(<3>; \l/) in L2(M, d ) for a given integer J 
naturally leads to a homogeneous wavelet system WS(\I/). Therefore, a homogeneous wavelet 
system WS(\I/) can be regarded as the limit of a sequence of nonhomogeneous wavelet systems 
WSj($; \&) as j — > — oo. Due to Theorem [7J we see that a nonhomogeneous orthonormal wavelet 
basis has natural connections to refmable structures. This motivates us to introduce the notion 
of a homogeneous wavelet system with the quasi-refinable structure. For a given homogeneous 
wavelet system WS(\I/) being a frame or a Riesz basis in L 2 (M. d ), we say that the homogeneous 
wavelet system WS(\I/) has the quasi-refinable structure, if there exists a subset $ of L 2 (M. d ) such 
that the nonhomogeneous wavelet system WSo($; is a frame or a Riesz basis in L 2 (M. d ) so that 
the homogeneous wavelet system WS(\I/) is its limit system. Lastly, we mention that results in 
this section for nonhomogeneous wavelet systems in L 2 (M. d ) can be generalized to more general 
function spaces such as Sobolev spaces and Besov spaces. For the case of Sobolev spaces, see [2TJ 
Theorem 7] and [21] for more detail. 

3. Frequency-based Nonstationary Wavelet Systems in the Distribution Space 

To characterize nonhomogeneous dual or tight wavelet frames in L2(M. d ), we shall take a frequency- 
based approach by studying frequency-based nonstationary wavelet systems in the distribution 
space. More precisely, we shall introduce and characterize a pair of frequency-based nonstationary 
dual wavelet frames in the distribution space. As pointed out in [21] for dimension one, such a 
notion allows us to completely separate the perfect reconstruction property of a wavelet system 
from its stability in various function spaces. Results in this section will serve as our basis to study 
nonhomogeneous and directional nonstationary tight wavelet frames in L2(M. d ) in the next section. 

Following the standard notation, we denote by 2t(M. d ) the linear space of all compactly supported 
C°° (test) functions with the usual topology, and ^'(M. d ) the linear space of all distributions, that 
is, 2t'(R d ) is the dual space of 3$(M, d ). By duality, it is easy to see that translation, dilation 
and modulation in (11. ip can be naturally extended to distributions in f§^(R ). For a tempered 
distribution /, by the definition of the notation fu,k,n m P-ip ? we have 

fu-,k,n = e~ ikn /(c/T)-i._ nik and fu- k = /(j/T)-i ; o,k- (3-1) 

In this paper, we shall use boldface letters to denote functions/distributions or sets of func- 
tions/distributions in the frequency domain. 

By L l ° c (M. d ) we denote the linear space of all measurable functions / such that J K \f(x)\ p dx < oo 
for every compact subset K of M. d with the usual modification for p = oo. Note that L l ° c (M. d ) is 
just the set of all locally integrable functions that can be globally identified as distributions, that 
is, L[ oc (R d ) C @'(R d ). For f G @(R d ) and if) e L[ oc (R d ), we shall use the following paring 

(f, </>>:=/ f(£MlK and (if,, f> := (f^0> = / ^(OWH- (3.2) 

When f G @(M. d ) and ip G @'(R d ), the duality pairings (f , rp) and (if), f) are understood similarly 

as (f , if,) := (^f) := W) ■ 

Let J be an integer and Nj, j J be d x d real- valued invertible matrices. Let $ and j ^ J 
be subsets of distributions. A frequency-based nonstationary wavelet system is defined to be 

oo 

FWS J (*;{* i }^ 5j ) = {¥> Nj; o,k : kGZ d ,<^G$}u|J{V>N j; o,k : kGZ^G*,}. (3.3) 

j=J 
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For the particular case Nj = l\P and SE^ = for all j ^ J, a frequency-based nonstationary wavelet 
system in (13.31) becomes a frequency-based nonhomogeneous (stationary) H-wavelet system: 

FWSj($; *) = {cp H j. 0jk : k G Z d , ip G $} U {^ ; o,k : j ^ J, k G Z d , x/> G (3.4) 

For a nonhomogeneous M-wavelet system WSj($; such that all the generators in $ and \I/ are 
tempered distributions, by ( 13. ip . the image of the nonhomogeneous M-wavelet system WSj($; \I / ) 
under the Fourier transform simply becomes the frequency-based nonhomogeneous (M T ) ^-wavelet 
system FWSj(<I>; where 

$ = {0 : g $}, * = : V e (3-5) 

For analysis of wavelets and framelets, as argued in [2T] for dimension one, it is often easier 
to work with frequency-based wavelet systems FWSj($; instead of space/time-based wavelet 
systems WSj($; \l/), though both are equivalent to each other under the framework of tempered 
distributions. Since we consider frequency-based wavelets and framelets in the distribution space 
@'(R d ), it is natural for us to consider FWSj(<£; *) C @'(R d ). 
Let Nj, j ' ^ J be d x al real- valued invertible matrices. Let 

# = {<pV--,V r } and $ = (3.6) 

and 

= . . . , ^} and = {ft' 1 , ft' 8 *} (3.7) 

be subsets of @'(R d ) for all integers j ^ J, where r, s G N U {0, +oo}. Let FWSj($; be 
defined in ( 13. 4 p and FWSj($; {^^fLj) be defined similarly. Generalizing the notion in (21] from 
dimension one to high dimensions, we say that the pair 

(FWSj(#; {*i}f =J ), FWS,(#; {*;}^)) (3.8) 

is a pazr of frequency- based nonstationary dual wavelet frames in the distribution space 3>'{R d ) if 
the following identity holds 

r oo s j 

EE( f '4o,k><4^)+EEE( f '<o,J«^g) = (2^) d (f, g ) (3.9) 

for all f , g G ^(R d ), where the infinite series in ( 13. 9 p converge in the following sense 

(i) for every f , g G S)(R d ), all the series 

r s j 

EE( f ^N j; o, k )(^N j; o, k ,g) and ^E( f '<;0,k)(^ ; o, k) g) (3-10) 

converge absolutely for all integers j ^ J; 

(ii) for every f , g G @(R d ), the following limit exists and 

^ m (E E< f -C>^*g> + EE E( f -0(*U-g>) = (2^) d (f,g)- (3.H) 

^ °° i=i wez d j=J t=\ wez d 

We say that the pair in ( 13. 8 p is a pair of frequency-based nonstationary dual wavelet frames in 
L2(M. d ) if (i) all elements in the two systems of the pair belong to L2(R ), (ii) each system in the 
pair is a frame in L2(R d ), and (iii) (13. 9p holds for all f , g G L2(R d ) with the series converging 
absolutely. It is straightforward to see that a space/time-based pair (WSj($; ty), WSj($; \I / )) is a 
pair of nonhomogeneous dual M-wavelet frames in L 2 (M. d ) if and only if the frequency-based pair 
(FWSj($; FWSj(<I>; 4f)) is a pair of frequency-based nonhomogeneous dual (M T ) ^-wavelet 
frames in L2(R d ), where are defined in ( 13. 5 p and 

$ = {0 : </> G * = {^:^G^}. (3.12) 
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Before we study and characterize a pair of frequency-based nonstationary dual wavelet frames 
in the distribution space @'(R d ), by the following result, we see that the above notion allows us 
to completely separate the perfect reconstruction property of a wavelet system from its stability 
in the function space L 2 (R d ). 

Theorem 9. Let Nj,j ^ J be d x d real-valued invertible matrices. Let <fr, <fr and ^fj,^fj be at 

most countable subsets of distributions on R d for all integers j ^ J. Then the pair in (13.81) is a 
pair of frequency-based nonstationary dual wavelet frames in L 2 (R d ) if and only if 

(i) there exists a positive constant C such that 

oo 

£^|(f^ Nj;0 , k >| 2 + E £ EK f ^N 3; o, k )| 2 ^^l|f|lL( K ,) ; Vf e®(R d ), (3-13) 
kez d j=J Ve*j kgz d 

oo 

Yl i(g^N J: o,k>r +E E E i(g^w>r < ci| g f i2(Rd) , v g e (3.w) 

(ii) i/ie pazr m (13.81) is a pair of frequency-based nonstationary dual wavelet frames in the 
distribution space ^'(R d ). 

Proof. The necessity part is trivial, since &{R d ) C L 2 (R d ). It suffices to prove the sufficiency 
part. By (13~T3|) . we have |(f, ^ N/;0 ,o)| 2 s$ C||f \\l 2m for all f G 3>(R d ). Therefore, we see that 



(•) Vn j; o,o) can be extended into a bounded linear functional on L 2 (IR d ), from which we conclude 
that Vnj;o,o can ^ e identified with a function in L 2 (R d ). Since Nj is invertible, we have ip G L 2 (IR d ) 
for all ip G 3>. By the same argument, we see that both systems in (13.81) have all their elements 
in L 2 (R d ). Consequently, since S3(R d ) is dense in L 2 (R d ) and both systems are countable sets, by 
a standard argument using density, we see that (13 . 13[) and (I3.14p hold for all f , g G L 2 (R d ). By 
item (ii), holds for f , g G @(R d ). Since ( ES} and (ETHj) hold for all f , g G L 2 (R d ), using 
Cauchy-Schwarz inequality, we see that (USD holds for all f , g G L 2 (R d ) with the series converging 
absolutely. Now by (13. 9p and using Cauchy-Schwarz inequality, it follows from (I3.13P and (13. 14j) 
that the left-hand sides of the inequalities in (13 . 13[) and (13. 14|) are no less than (27r) d C -1 ||f || 2 2 ( R d) 
and (27r) d C _1 ||g||^ 2 , IR( i- ) , respectively. This completes the proof of the sufficiency part. □ 



A pair of frequency-based (or space/time-based) nonstationary dual wavelet frames can be 
generalized from the function space L 2 (R d ) to a pair of dual function spaces (SB, SB') where SB is 
a Banach function space contained in ^'(IR d ) and SB' is its dual space. For example, S3 can be a 
Sobolev or Besov space. The result in Theorem [9] holds for such more general function spaces by 
replacing item (i) with suitable stability condition, more precisely, the boundedness of the primal 
system in the function space SB and the boundedness of the dual system in the function space SB' . 
That many classical function spaces can often be characterized by the wavelet coefficient sequences 
largely lies in the fact that we have the perfect reconstruction property in item (ii) while the two 
wavelet systems are often bounded in many function spaces. For the case of Sobolev spaces, see 
[21| Theorem 7] and |24] for more details on pairs of dual wavelet frames in function spaces other 
than L 2 (R d ). 

To characterize a pair of frequency-based nonstationary dual wavelet frames in the distribution 
space, let us first look at the absolute convergence of the series in (I3.10p . By the following result 
which generalizes [21], Lemma 3] , we see that the absolute convergence of the series in (13 . 1 0[) holds 
under a very mild condition. 

Lemma 10. Let U be a d x d real-valued invertible matrix and if), if) G L l 2 oc (R d ). Then for all 

f,geS>(R d ), 

^(f,^ ; o,k)(^;0,k,g) = (2n) d I f (0g(£ + 2nU-'k) WO^t + 2vrk) d£ (3.15) 

kez d jRd kez d 
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with the series on the left-hand side converging absolutely. Note that the infinite sum on the 
right-hand side of (13.151) is in fact finite. 

Proof. Define two 27rZ d -periodic functions h and h as follows: 



h(0== ^f(f/- 1 (^ + 27rn))V(e + 2vrn) and h(£) := £ g^ 1 ^ + 2nn))^ + 2nn). (3.16) 

Since f,g G S$(M. d ), the infinite sums in ( 13 . 16f) are in fact finite for £ G [— 7r, ir] d . Moreover, by 
ip,ip E L l £ c (M. d ), we deduce that h,h 6 L 2 (T d ). By calculation, we have 

/ h(Oe**de = | det U\V 2 (f, ^ ; o,k) and / h(£K k ^ = I det f/| 1/2 (g, $ U;0 , k ). 

Now by Parsevel identity for periodic functions in L2(T d ), one can easily deduce that the left-hand 
side of (I3.15P is equal to 

(2n) d \detU\- 1 [ h(OW)d^ = (27r) d \detU\- 1 [ iiU^oWHO^ 

J{-n,Tr) d JM. d 

= (2n) d [ fiowmWln, 

JR d 

from which we see that (13.151) holds. Since f , g are compactly supported and U is invertible, 
f (f)g(f + 2ttU~ 1 V) = for all (Gl d provided that ||k|| is large enough. Hence, the infinite sum 
on the right-hand side of (13.151) is in fact finite. □ 

The assumption of membership in L^lW 1 ) is only used in this paper to guarantee (13.151) with 
the series on the left-hand side of (13.151) converging absolutely. 

The following result characterizes a pair of frequency-based nonstationary dual wavelet frames 
in the distribution space. 

Theorem 11. Let J be an integer. Let Nj, j ^ J be d x d real-valued invertible matrices such that 
A := U c *L J [N~ l Z d ] has no accumulation point and lim ||N,-k||=0 V k G A. (3.17) 

Let $ in (13. 6p and Vt^,^ in (13.71) be finite subsets of L l £ c {M) for all integers j #J J . Then 
the pair in (13.81) is a pair of frequency-based nonstationary dual wavelet frames in the distribution 
space 3!'(M. d ) if and only if 

j'-i 

and 

oo 

2j jk (NjO + E Z *f ( N ^) = °> a - e - ^ M d , k G A\{0}, (3.19) 

3=J 

(the infinite sums in ( 13. 19[) are in fact finite.) where 

r 

2* (£) : = ^Ipiig)? 1 ^ + 2nk), k G Z d and X% := k G R d \Z d , (3.20) 
i=i 

S 3 

2* .(£) := ^iftMfcWte + 2rrk), k G Z d and 1%. := k G R d \Z d . (3.21) 



lim (2g(N J -) + E J * ,(N 3 -),h) = (l,h) Vhe^(l d ) (3.18) 

3=J 
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Proof. For f , g G @(R d ), denote 

#(f,g) := £ ^(f,^ Nj;0 , k )(^N j; o, k) g) + EEE< f '<;o.k><<;0,k.g>- ( 3 -22) 

By Lemma fTUl we have 

#(f,g) = (27r) d /* ^f(Og(^ + 2 7 rk)(x^ k (N J + ^Xj/IN.O)^. (3.23) 

On the one hand, since f and g are compactly supported, there exists a positive constant c such 
that f(f)g(f + 2vrk) = for all f G M d and |k| ^ c. On the other hand, by our assumption in 
f)3.17p . {keA : |k| < c} is a finite set and therefore, there exists a positive integer J" such that 
Njk G" Z d for all j ^ J" and k G A\{0} with |k| < c. That is, under the assumption in ( 13. 17ft . for 
all J' ^ J", (13T23D becomes 

j'-i 



(3.24) 



#(f,g) =(2tt)<* [ f(0g(0(2S(Nj0 + X)^(NiO)de 

„ oo 

+ (2n) d / £ f(0g(e + 2vrk)(j^ k (N J + E^f (N^))^ 

jRd keA\{0} i=j 

If (13TT8D holds, then we deduce from (l3T24j) that 

f J'-i 

#(f,g) = (2n) d / f(Og(0(^(NjO + (N ; 0k- (3-25) 

Now it follows from ( 13.181) that 



lim Sj(f,g) = (27r) d (f,g>, f,gG^(R rf ). (3.26) 

J'— !-+oo 

Therefore, the pair in (I3.8P is a pair of frequency-based nonstationary dual wavelet frames in the 
distribution space @'(R d ). 

Conversely, ( I3.26P holds. By our assumption in (I3.17p . (I3.24p holds for all J' ^ J". Note that 
the set A is discrete and closed, that is, for any k G A, e k := infygA.Uk} \\u ~ Ml/2 > 0. For any 
k G A\{0} and £ e we deduce from ( l3T24"j) and ( 13T26|) that for J' ^ J", 

» oo 

(2n) d / f(0g(e + 2vrk)(x^ k (N, / + Vl^N^))^ = Sf (f,g) = (2n) d (f,g) = (3.27) 

for all f, g G ^(M d ) such that suppf C B £y (^ ) and suppg C B £k (£ - 27rk), where -B £k (£o) : = 
{£ G M d : ||£ — £ || ^ ^fc} is the ball with center £ and radius e^. Now we can easily deduce from 
the above relation in (I3.27P that identity in (I3.19P holds for almost every £ G B ek (£ ). Since £ 
is arbitrary, we conclude that (13.191) holds. By (13.191) and (13.241) . we see that (I3.25P holds, from 
which we see that (I3.18P holds. □ 

As we shall see in the proof of Corollary [151 the condition in (13.171) is automatically satisfied if 
Nj = l\P for all j ^ J and all the eigenvalues of N are less than one in modulus. Using a more 
technical argument, we also see that Theorem [11] (as well as other results in this paper) still holds 
if the assumption, that 3>, V&j, j are finite subsets of L l 2 C (W d ) for all j ' ^ J in Theorem [TTl is 
replaced, for example, by the assumption 

E m-)\ 2 , E m-)\ 2 , E i-0(-)i 2 , E i-0(-)i 2 , j > j 

VG* VG*j */>S*j 
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belong to L l 2 oc (R d ). 

The following result is a direct consequence of Theorems [9] and [TT1 

Corollary 12. Let Nj, j ^ J be d x d real-valued invertible matrices. Let $ and ^/j be subsets of 
^'(IR d ) for all j ^ J. Then the following statements are equivalent to each other: 

(i) FWSj(<fr; is a frequency-based nonstationary tight wavelet frame in L 2 (R d ), that 
is, *j C L 2 {R d ) for all j ^ J and 

oo 

£^|<f,<p Nj; o, k >| 2 + E £ EK f ^N 3; o, k )| 2 = (2vr) d ||f||i 2(Rd) , VfGL 2 (f d ); (3.28) 
ye* kez d i=J kgz d 

(ii) (FWSj($; {^-^j), FWSj($; {^'j^j)) is a pair of frequency-based nonstationary dual 
wavelet frames in the distribution space @'(R d ). 

If in addition ( 13.17!) holds and all j ^ J are finite subsets of L l 2 oc (M d ) , then any of the above 

statements is also equivalent to 

(hi) (13.181) and (13.191) are satisfied with $ = $ and = ^ for all integers j ^ J. 

Proof. By TheoremEl (i)=^(ii). To see (h)^(i), since (FWSj($; FWSj($; {*j}£Lj)) is 

a pair of frequency-based nonstationary dual wavelet frames in the distribution space, by definition, 
we can easily deduce that (13.131) must hold with C = (2n) d . Now by Theorem [9j (ii)=>(i). The 
equivalence between items (ii) and (iii) follows directly from Theorem [TTJ □ 

We have the following result on a sequence of pairs of frequency-based nonstationary dual 
wavelet frames in the distribution space. 

Theorem 13. Let Jq be an integer. Let Nj,j ^ Jo be d x d real-valued invertible matrices. Let 

<l>, = W h \ ■■■■ <P*' ri h <I> ; = {&'\ •••• ^ M 1 (3.29) 
and ^fj, in (13. 7j) be subsets of 2)'(Jsi d ). Then the following statements are equivalent: 

(i) the pair (FWSj($j; j), FWS j(l j; j)) is a pair of frequency-based nonsta- 
tionary dual wavelet frames in the distribution space £^'(IR d ) for every integer J ^ Jo; 

(ii) the following identities hold: 

lh ? E E( f '^ ; o, k )(^;0. k) g) = (2vr) d (f,g), Vf,g G (3-30) 

'I— >+oo < * ■'— ' J ' ' ' 

«=i kez d 

and 

E E < f > <o.k> (^,k. + E E ( f > <;o,k> <^ ; o, k > g) 
€=1 k ez d 

witt a// i/ie series converging absolutely. 

If in addition (13.171) holds and all ^ Jo are finite subsets of L l 2 oc (R d ), then any 

of the above statements is also equivalent to 

(iii) the following identities hold: 
.lim (^^(N,-)^(N,.),h\ = (l,h), Vhef(l d ) (3.32) 



-oo 

1=1 



16 BIN HAN 

and for all integers j ^ Jo, 

^f(NiO-4f(N,0 = ^(N^iO, a.e. £ G M d , k G [Nj 1 ^ U [N7 + \Z d ], (3.33) 
where ,k 6 l d are defined in (I3.2ip and 
r j 

4 J (0:=E^)^ + 27rk )' kGZ " anrf J *, :=0 kGK d \Z d . (3.34) 

£=1 

Proof. (i)=^(ii). Considering the difference between two pairs at consecutive levels j and j + 1, 
from (13.111) we see that (13.311) holds. Now by (13.311) . it is straightforward to see that 

rj J'-l s j 

= EE( f ^N;f ; o >k )(^ ; o 1 k>g}> > 

Therefore, (T3T3TH) holds. 

(ii)=^(i). Since all the series in (I3.3ip converge absolutely, all the series in (I3.10p converge 
absolutely. To see item (i), by (13.311) . we deduce that ( I3.35|) holds. Now (13.111) with <3? = and 
$ = $j follows directly from (13T30D and fl3T35|) . 

(ii)<^(iii). By Lemma [TUl we see that (I3.3ip is equivalent to 

f £ f (£)g(£ + 27rk) (Z^f (N.O + 4f (N,0 - 2j[ j; k (N i+1 0) # = 0, (3.36) 

where Aj := [N~ 1 Z d ] U [NjV^Z^]. Since Aj is discrete, by the same argument as in Theorem [TTj we 
see that (13.361) is equivalent to (I3.33[) . 

By Lemma [TDJ, we see that ( 13.30[) is equivalent to 

lim / Yl f (Og(C + 27rk)X^ k (N,e)rfe = <f, g>, f.ge*. (3.37) 

By our assumption in (13.171) . as in the proof of Theorem [TTJ there exists a positive integer J" such 
that f(f)g(f + 27rk)i^ k (NjO = for all £ G IR d , k G A\{0} and j ^ J" . Consequently, for j ^ J" , 
( 13.371) becomes 

.lim / f(0i(02» J (N i 0de= / f(0i(0de, f,g G^(M d ), 



which is equivalent to (I3.32p . □ 
We point out that (13.331) can be equivalently rewritten as 

4,(0 +2^(0 =2j^ N7lk (N i+1 N i - 1 0, a.e. ^l d ,ke Z d U [(N^Nt 1 )-^]. (3.38) 

Due to its similarity to MRA (multiresolution analysis) structure ([L2J ETJ EHj), the relations in 
(13.31 j) (or more general, (I3.35P ) could be called the quasi-MRA structure. One important property 
of dual wavelet frames is its frame approximation order (see [H]). Due to ( 13.351) . we observe that 
the frame approximation order associated with the pair in Theorem [13] is completely determined 
by the generators <E»j, 4fj,j ^ Jo and it has nothing to do with ^ j, *S?j,j ^ Jo- 
The following simple result on matrices will be needed later. 



NONHOMOGENEOUS WAVELET SYSTEMS IN HIGH DIMENSIONS 



17 



Lemma 14. Let M be a d x d real-valued matrix. Assume that Ai, . . . , A<f are all the eigenvalues 
of M ordered in such a way that |Ai| ^ |A 2 | ^ ■ ■ • ^ |A d |. Define A min := |AjJ and A max := |Ad|. 
Then for any e > 0, there is a norm || ■ ||m on ~R d such that 

(X min -e)\\x\\ M ^\\Mx\\ M ^(X max + e)\\x\\ M , VigM". (3.39) 

Moreover, if every eigenvalue X of M satisfying |A| = A min /ias t/ie same algebraic and geometric 
multiplicity, then A min — e on the left-hand side of (13.391) can be replaced by A min . Similarly, if 
every eigenvalue X of M satisfying |A| = A max nas the same algebraic and geometric multiplicity, 
then A max — e on the right-hand side of (13.39!) can be replaced by A max . In particular, if M is 
isotropic (that is, M is similar to a diagonal matrix with all the diagonal entries having the same 
modulus), then ||Mx||m = | det M| 1//d ||a;||M for all x £ IR d . (Conversely, this identity implies that M 
is isotropic.) 

Proof. There exists a d x d complex-valued invertible matrix E such that EME" 1 is the Jordan 
canonical form of M. Denote 

ro i ••• o" 

o o ••■ 1 
p o ••• o_ 

Let XI r + F be a Jordan block in the Jordan canonical form EME^ 1 of M. For any e ^ 0, we have 
diag(l, e -1 , e~ 2 , . . . , e 1 " r )(A/ r . + -F)diag(l, e, e 2 , . . . ,e r ~ 1 ) = XI r + eF. Therefore, XI r + F is similar 
to XI r + eF. Consequently, there exists a d x d complex-valued matrix G such that GMG' 1 is a 
diagonal block matrix with each block being the form of XI + eF. Define a norm || • ||m on C d by 

—T 

\\ x \\w\ := ll^ x l| 2 := % T G Gx. Restricting the norm || • ||m on M d , now we can easily deduce that 
(ET351) holds. □ 

As a direct consequence of Theorem [121 we have the following result on a pair of frequency-based 
nonhomogeneous dual wavelet frames in the distribution space. 

Corollary 15. Let M be a d x d real-valued expansive matrix and define N := (M T ) _1 . Let <&, $ 
in (I3.6P and 

* = {V 1 ,...,^"} and * = {i/? 1 ,...,^} (3.40) 
be finite subsets of L l 2 C (¥L d ). Then the following statements are equivalent: 

(i) (FWSj(<fr; St), FWSj(<fr; St)) is a pair of frequency-based nonhomogeneous dual H-wavelet 
frames in the distribution space 3t'(M. d ) for some integer J; 

(ii) item (i) holds for all integers J; 

(iii) for all f , g £ £t(M. d ), the following identities hold: 

r 

lh ? E E( f '^;0, k )(^;0, k ,g) = (2vr) d (f,g> (3.41) 



and 



(3.42) 



EE^' ^N;0,k)(^N;0,k>g); 

t=l kez d 



(iv) lim^+oo Y^e=i ¥^(N J '-)y (I\P-) = 1 in t/ie sense of distributions and 

4(0 + 4(0 = (NO, a.e. f £ M d , k £ Z d U [N -1 Z <i ], (3.43) 
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where l\ is defined in ( 13 . 20 j) and 



X£ (f ) : = V>'(0^(£ + 2vrk), k G Z d and X% := 0, k G M d \Z d . (3.44) 



Proof. We first show that the condition in (13.171) is satisfied. Since Nj = l\P and all the eigenvalues 
of N are less than one in modulus, by Lemma 03} it is trivial to conclude that lim^oo ||N J x|| = 
for all x G M. d . By Lemma [T41 again, it is also trivial to see that the set A := U^Lj[N~ 3 °Z <i ] = 
U£L < /[(M T ) J 'Z <i ] is discrete. So, the condition in (13. 1 7[) is satisfied. 

Now all the claims in Corollary [TS] follow directly from Theorem [I3J □ 

As another application of the notion of a pair of frequency-based nonstationary dual wavelet 
frames in the distribution space, we have the following result which naturally connects a nonsta- 
tionary filter bank with wavelets and framelets in the distribution space. 

Theorem 16. Let {Mj}°? =1 be a sequence of d x d integer invertible matrices. Define N := Id 
and Nj := (Mf • • • Mj") -1 for all j G N. Assume that flBTTTj) holds. Let a, and a,-, j G N be 2irZ d - 
periodic measurable functions such that there exist positive numbers e,e,Cj,Cj,j G N satisfying 

|1 - a,(N,OI < CM\\ e , |1 " a,(N,OI < Q||£|| £ ~, a.e. £ G R d ,j G N (3.45) 

and 

oo oo 

y~]Cj <oo and y~]Cj < oo. (3.46) 

i=i j=i 

For all j G NU {0} ; define 

oo oo 

<Pj{® ■= II ^(Ni+nNj^) and := ^(N^Nj^), £ e (3.47) 

n=l n=l 

Taen a// <Pj,j G NU {0} are well-defined functions in L l ^(M. d ) satisfying 

^-i(Mje) = a,(e)^(0 and ^ i _ 1 (MjO = a,-(0^(0> «.(eR d ,jeN. (3.48) 

Let ij:L , . . . , b A1 , . . . , b ijSj _ 1 and A1 , . . . , 6^ r ._ l} b jtl , . . . , b,, Sj _ i; un't/i r^i, G N U {0} 

and j G N, be 2it'L d -periodic measurable functions in L l 2 C (M. d ). For j G N, define 

ift-^ii) := 0,AO^ dO and ^{£) := * = 1, . . . , r,-_ lf (3.49) 

^'-^(MJO := b^(0Vi(0 and ^(MjO := b^(0^'(0, € = 1, . . . , (3.50) 

Tnen m f l3T2TJ|) ana 1 m f l3"Tj) are stz&sete o/L z 2 oc (M) /or a// j G NU {0}. Moreover, 

the pair 

(FWSj(# j; {^}f =J ), FWSj(* j; {^}f =J )) (3.51) 

is a pair of frequency-based nonstationary dual wavelet frames in the distribution space @'(M.) for 
every integer J ^ 0, if and only if, for all j G N, 

s j-i 

e 3 {Mj^A0^(0 + ^>^(£)M0 = 6 i+1 (0, a.e.£ G (7 w n^ (3.52) 

Sj_l 

e^MjO^COa,^ + 2tto;) • )!,,•,(£ + 2™) = 0, a.e. £ G o>. n (^. - 2™) (3.53) 

/or a// w G fi Mj \{°} Wlth ■= [(Mj)- 1 Z d ] n [0, l) d ; and 

lim 0j + i(Nj-) = 1 in the sense of distributions, (3.54) 

j— >+oo 
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where 



kez d kez d 



rj-i 



a Vi := {?Gl d : \ipj(Z + 2vrk)| + o}, := K d : ^ + 2vrk)| ^ o} (3.55) 

and 

£=1 

Proof. Using the inequality \z\ ^ e' 1-2 ' for all z G C, by (I3.45p . we deduce that 

m 2 

| n a J+ „(N, +n 0| < I^^OW)! <; e ll*ll'£3»*GH*. ^ e o|| € ||- (3 57) 

n=mi 

for all m!,m 2 G N with mi < m 2 , where C := J^Li^j < 00 • O n the other hand, we have the 
following identity 

m2 m2 ni2 

!- II a i+"( N i+«0 = XI (i - a i+m( N i+mO) ( a J+n( N J+"0) > 

n=mi m=mi n=m+l 

where we used the convention n™=m 2 +i := ^ Therefore, by f!3.45j) and (I3.57p . we have 

m 2 m 2 m 2 

1 - J] a, + „(N, +n 0| < e c ll«" £ 2 |1 - a J+m (N, +m OI < ^"1^ £ C i+m . (3.58) 

n=mi m=mi m=mi 

Since C = Q < oo, the above inequality implies the convergence in L^ c (IR d ) of the infinite 

product n^Li a i+"( N i+"0- Since ¥>j(NjO = Iir=i a i+n(Nj+nO> it follows from (13.571) that ipj G 
L^ c (R d ). Since all 0^,b^ G L 2 oc (R d ), it is evident that all tp?*,^* are elements in L 2 oc (R d ). 
Similarly, we can prove that (f>j G L^f(IR d ) and all ipi^,^^ are elements in L 2 oc (IR d ). 

By (I3.47p . we have y?j_i(Nj_i^) = aj(Nj-£)<^(l\l 3 -£), from which we can easily deduce that (13.481) 
holds. Since MjN^ = Nj_i and M., is an integer matrix, we see that (I3.33j) with j being replaced 
by j — 1 is equivalent to 

r j-i s j-i 

^- 1 ^(MjO^" 1/ (Mj(e + 27rk + 2nco)) + ^'-^(Mjf )^'~ M (Mj(£ + 2vrk + 2ntu)) 

1=1 rj (3.59) 

= 5(u)J2^p^)(p j ' e ^ + 2vrk + 2ttu) 

for almost every £ G M d and for all k G Z d and u; G fluj, where 5(0) = 1 and 5(oj) = for all 
w^O. Now by (I3.49p and (I3.50p . it is easy to deduce that (I3.59P is equivalent to (13. 52j) and (I3.53p . 
By (I3.58p . we deduce that 

oo oo 

\1 - <p,m)\ ^ e™^Z\\ e T, C m and |1 - ^(N^)| < e^'llflf £ C m , (3.60) 

m=j+l m=j+l 

where C := Y^jLi^j < 00 • By (I3.46p . the above inequalities imply that for any bounded set K, 
there exists a positive integer J% such that 

\ < IVi(NiOI ^ | 5 < l^(NiOI < ~, a.c e G X, j ^ Jjf. (3.61) 
It follows from fl£gg} that 

X^(N ; 0^!N ; = i+1 (N i O^(N5O'^(N i O- 
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By (13. 60 p . (I3.6ip and the above identity, using Lebesgue dominated convergence theorem, we 
conclude that lim^oo YltLi P 3 ' £ (^j')'P : '' i {^j') = 1 m the sense of distributions if and only if (I3.54p 
holds. The proof is now completed by Theorem [13] □ 

The conditions in (I3.52p and (I3.53P in Theorem [TBI generalize the well-known Oblique Extension 
Principle (OEP) in [5j H31 E] (also see [151 E3 EQ] [22]) from homogeneous wavelet systems to 
nonhomogeneous and nonstationary wavelet systems. The stationary case of OEP with rj = 1 
and Oj i = Oj^i = 1 is given in [29j [30] for homogeneous wavelet systems. Let M be a d x d integer 
expansive matrix and N := (M T ) _1 . In Theorem [T6] if Mj = M, a.j = a, a,- = a for all j G N and if 
there exist positive numbers e, e, C, C such that 

|l-a(OKC||e|| e and |1 - a(0| < C||£f, a.e.£ eR d , (3.62) 

by Lemma HU then it is easy to see that (I3.17p . (I3.45p . and (I3.46P are satisfied. 

4. Nonhomogeneous and Directional Nonstationary Tight Framelets in L 2 (M. d ) 

In this section, we shall study and construct nonhomogeneous (stationary) tight wavelet frames 
with a minimum number of generators by proving Theorem [2j By modifying the construction in 
Theorem [2] we also show that we can easily construct directional nonstationary (more precisely, 
nonhomogeneous quasi-stationary) tight wavelet frames in L 2 (M d ). Moreover, both constructed 
systems are associated with filter banks and can be derived by Theorem [16] with such filter banks 
if the underlying dilation matrix is an integer expansive matrix. 

As a direct consequence of Corollary [12] and Theorem [TBI we have the following result. 

Corollary 17. Let M be a d x d real-valued expansive matrix and define N := (M T ) _1 . Let J 
be an integer. Let <fr in (13. 6p and SE^ in (13. 7p be finite subsets of L l 2 oc (M. d ) for all j ^ J . Then 
FWSj(<fr; is a frequency-based nonstationary tight wavelet frame in L2(M. d ) for every 

integer J ^ Jq, that is, \&j C L2(M. d ) for all j ^ Jq and 

r oo Sj 

E E K^WI 2 + EEK f 'C,J| 2 = (2vr)1f||i 2(Rd) , Vf G L 2 (M. d ), J ^ J , (4.1) 
e=i kgz d j=J e=i 

if and only if, 

r 

lim \ \cp e (W-)\ 2 = l in the sense of distributions (4-2) 
l=\ 

and for all j ^ J 0; 

r Sj 

E^V(£ + 2vrk) + E^W^(£ + 2vrk) 

e=i 

r 

E^W)</(N(£ + 27rk)), a.e.£ G R d ,ke Z d n [N _1 Z d ], 



e=i i=i 

r 



(4.3) 



E ^(0^(C + 2vrk) + E ^ l ^W'\i + 2vrk) = 0, a.e. £ G M d , k G Z rf \[N" 1 Z' i ], (4.4) 

£=1 f=l 
r 

E ^(NO^(N(e + 2vrk)) = 0, a.e. £ G R d , k G [N _1 Z d ]\Z d . (4.5) 



Moreover, if the following property holds: 



h(Oh(f + 2?rk) = a.e. £ G M d , k G Z d \{0}, he$U (U^ Jo ^), (4.6) 
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then all the conditions in CI4.3[) . (I4.4p . and (14.51) are reduced to the following simple condition 

r s j r 

E l*>'(OI 2 + E l^(0l 2 = E l^( N OI 2 , o.e. £ G M d , J ^ Jo- (4.7) 

1=1 1=1 i=\ 

Proof. By the proof of Corollary [TS] we see that the condition in (13.1 7\i is satisfied. It is trivial to 
check that (13.331) with := <&j := <E» and SE^ := is equivalent to all the conditions in (I4.3p . 
(I4.4p , and (14.51) . Note that the assumption in (14.61) simply means X^, (0 = and 1\ (0 = for 
almost every £ G K d , k G M d \{0}, and j ^ J . Now it is easy to see that all the conditions in 
(14. 3p . (14.41) . and (14. 5 p are reduced to the simple condition in (14. 7\\ . Since $>j = = and rj = r, 
we have Y2=i <P'-M N (N,-£) = ££ =1 |</(l\P'0| 2 - Now it is straightforward to see that fl3T32|) 
is equivalent to (14. 2 p . The proof is completed by Theorem [13] and Corollary [12] □ 

The sequence of systems FWSj(<E>; {SE^-j-^j), J ^ J in Corollary [17] is of particular interest 
in applications. The only difference between these systems in Corollary [T7J and nonhomogeneous 
wavelet systems lies in that the sets ^Sfj of wavelet generators depend on the scale level j and 
therefore, the wavelet generators are not stationary. This freedom of using different sets SE^ at 
different scales allows us to have different ways of time-frequency partitions at different resolution 
scales and therefore such a freedom is of importance in certain applications. Since these systems in 
Corollary [T7J are between nonhomogeneous stationary wavelet systems and (fully) nonstationary 
wavelet systems, we shall call them nonhomogeneous quasi-stationary wavelet systems. 

Now we are ready to prove Theorem [2] in Section 1. 

Proof of Theorem^ Let A m i n be the absolute value of the smallest eigenvalue of M in modulus. 
Note that A m i n > 1. Define B c (0) := {£ G IR d : ||0Im t < c } ; where || ■ ||m t is the norm in Lemma HH 
with M being replaced by M T and with an e > satisfying A := A m i n — e > 1. By (I3.39p . we see 
that B c (0) C B Xc (0) C M T B c (0). Choose A such that < A < 1. For small enough c, we have 

B c (0) C B Xc (0) C M T B c (0) C (-A tt, \ n) d and (M T ) 2 B c (0) C (-A vr, A 7r) d (4.8) 

For any p > 0, there is an even function h p in C°°(M) such that h p takes value one on (— p/2, p/2), 
is positive on (— p, p), and vanishes on M\(—p,p). Such a function h p can be easily constructed, 
for example, see rTS] Section 4]. For sufficiently small p > 0, we define 

. XB {1+x)c/2 ( )*h p (\\ ■ ||) 

For sufficiently small p > 0, it is evident that is an even function in C°°(M. d ), ^ ¥?(0 ^ 1 for 
all £ G M d , and 

<p(£) = 1 V £ G 5 C (0) and y>(0 = V £ G R d \(M T 5 c (0)). (4.10) 

Now we define 

i>(0--=VW) with 77(0 := MN£)) 2 - MO) 2 , eet d , (4.n) 

where N := (M 7 ; -1 . By (14. 8 p and (14.101) . we see that r/ vanishes on B c (0) and is supported inside 
(M T ) 2 J B c (0)C(-A 7r,A vr) d . 

Since (p G C°°(R d ), it is obvious that 77 G C°°(R d ). Since < cp < 1 and by Olll . 

V?(N0 = 1 whenever <p(£) ^ 0. (4.12) 
We deduce from (I4.12p that ^ r)(£) ^ (cp(NO) 2 - Therefore, ij) is well defined and 

< V>(M T < <P(Z) < h and V(M T = V(0, Z G K d \£ c (0). (4.13) 

To prove that i/> G C°°(R d ), we first show that if v?(£o) = or 1, then all the derivatives of <p at 
£o must vanish. In fact, if <p(£o) = 0, then by the definition of <p in (14 .9p we see that the supports 
of h p (\\ ■ ||) and Xb, 1+X)c/2 (o)(£q — ■) are essentially disjoint. Consequently, it follows from (14. 9p that 
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all the derivatives of (p must vanish at £ - Similarly, if <^(£o) = 1> then by (14. 9p we see that the 
support of h p (\\ ■ ||) is essentially contained inside the support of XB( 1+a)c / 2 (o)(£o — ■)• Consequently, 
it follows from (14. 9[) that all the derivatives of cp must vanish at £o- 

Now we show that if) £ C°°(IR d ). In fact, if ry(£o) 7^ 0, since 77 G C°°, then it is trivial to see 
that if) = ^Jfj is infinitely differentiable at £o- If *7(£o) — 0, by (14. lip and (I4.12p . we must have 
either (p(£o) = y(N£o) = or <p(£o) = ^(N^o) = 1- By what has been proved, we see that all the 
derivatives of 77 must vanish at £o- Using the Taylor expansion of 77 in a neighborhood of £0, now 
it is easy to see that if) = ^/fj must be infinitely differentiable at £ with all its derivatives at £ 
being zero. Therefore, if) G C°°(lR d ). 

Hence, we constructed two functions ip,if) <E C°°(M. d ) vanishing outside [— n,7i] d and 

I^(0| 2 + I^(0| 2 =IV(N0| 2 V(Gt d . (4.14) 

Define and if) via <fi := cp and if) := if). By Corollary [T7] all items (i)-(iii) of Theorem [2] hold. 
Now we define 27rZ d -periodic functions a k , k G Z d by 

a k (0 := e- lk - MT ^(M T 0, £ G [-vr, n) d , k G Z d . (4.15) 

By dHl2]), (p(£) = 1 whenever v?(M T £) ^ 0. Since cp(M T -) vanishes outside B c (0) C (-A tt, \ n) d , 
now it is straightforward to see that a k G C°°(T a! ) and a k (£)v?(£) = e _ik-MT V(M T for all £ G M d . 

Let h Ao be a C°° function such that h\ = 1 on [— Aovr, Aovr] d and h\ = on [— n, — 7r) d \[— (Aq + 
Eo)ir, (A + £o)^] d , where < e < 1 — A . Such a function h Ao can be easily constructed using 
tensor product (see [TU Section 4]). Now we define 27rZ d -periodic functions b k , k G Z d by 

. = fe- kMT ^ (0^, ^ and £ G [-*,*)'; 

\e- kMT «h Ao (0, ^(O = 0andeG[-7r,7r) d . 

By fT4TT3D . it is trivial to see that e~ ik ' MT ^(M T £) = b k (£)v(0 for all £ G M d . Now we show that 
b k G C°°(T d ). Note that b k is infinitely differentiable at £ for all £ G [-vr, 7r) d such that ip(£) ^ 0. 
By the fact that h Ao vanishes outside [— (A + Sq)^, (A + £o) 1T ] d an d A + £0 < 1, by (I4.13p . we 
conclude that b k G C°°(T d ). Therefore, item (iv) holds. 

Note that (I1.7P is a direct consequence of Proposition |4j This completes the proof. □ 

When M in Theorem [2] is a d x d integer expansive matrix, the nonhomogeneous tight wavelet 
frames in Theorem [2] can be obtained via Theorem [161 from its filer bank (a ,b ) with 

Tj = 1, Sj = 1 and a.j = a,- = a , b^i = hj t i = b , = 0j \ = 1, j G Z. 

One potential shortcoming of the construction in Theorem [2] is that the support of if) has a ring 
structure and if) behaves like a radial basis function. In fact, the functions and if) can be con- 
structed to be radial basis functions for special dilation matrices such as M = 2J^, see [IS] for more 
details. However, as pointed out in [3J, for high dimensional problems, directionality of a transform 
is an important feature. All stationary wavelets and framelets have the isotropic structure and 
cannot capture very well singularities other than the point singularity in high dimensions. 

In the rest of this section, going from frequency-based nonhomogeneous tight wavelet frames to 
frequency-based nonhomogeneous quasi-stationary tight wavelet frames, we shall see that direc- 
tionality can be easily achieved by modifying the construction in Theorem [2] and using a simple 
splitting technique in [18]. The key lies in that ip has a very small support and if) is linked to <p 
via the relation in (I4.14p . By (14. 7p in Corollary [T7J the single generator if) in (I4.14p at every scale 
level can be easily split into many generators if)-'' 1 , ■ ■ ■ , if)^' Sj with directionality so that 

X> J 'W = l</>(0| 2 , a.e.£GM d . (4.17) 
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In the following, we describe the idea of the splitting technique by increasing angle resolution. 
For each integer j, pick a positive integer Sj and construct C°° functions (3j^, . . . , (3j )S . on the unit 
sphere of ~R d such that 



^|/3^(0| 2 = 1 VCGR d with ||e|| = l. (4.18) 



e=i 

To capture various types of singularities, preferably the significant energy part of each con- 
centrates near a point or an n-dimensional manifold on the unit sphere with ^ n < d. Let ip 
and i\) be constructed as in the proof of Theorem [21 Now, we can split ip by defining 

^(0:=^(0ft4), e=l,..., 8j and j G Z. (4.19) 

Since i\) vanishes in a neighborhood of the origin, all ip j > e G C°°(M a! ) and fl4T7D is obviously true. 
Consequently, for all j G Z, 

l^(0| 2 + £ |^(e)| 2 = l^(NO| 2 , £ e M d . (4.20) 

£=1 

Define *j := • • .,ip j ' Sj }. Now by Corollary Q7] and Theorem |2l FWSj({c^}; j) is a 

frequency-based nonstationary tight wavelet frame in L2(R d ) for all integers J. 

We mention that the main idea of the construction in Theorem [2] and the above splitting 
technique already appeared in [TS], Proposition 3.8 and Section 4] for M = 21^ and for homogeneous 
tight wavelet frames in L 2 (M. d ). 

Many directional nonhomogeneous quasi-stationary (and nonstationary) tight wavelet frames 
in L 2 (M <i ) can be constructed. To illustrate the above general procedure, using polar coordinates, 
we present the following result for the special case d = 2 and M = 2I 2 . In the following result and 
its proof, note that r and 9 are used for polar coordinates and a complex number xe %e is identified 
with the point (rcos#, rsin#), that is, M 2 is identified with the complex plane C. 

Theorem 18. Let m be a positive integer and ^ p < 1. Then there exist two real-valued 
functions <p and r\ in L2(M d ) fl C°°(IR 2 ) satisfying all the following properties: 

(i) (p is a compactly supported radial basis function in C°°(R 2 ) and there exists a 2^7? -periodic 
trigonometric function a in C°°(T 2 ) such that 0(2£) = a(£)0(£) for all £ G M 2 ; 

(ii) r] has the tensor-product structure in polar coordinates and there exist positive real numbers 
r i) r 2, 9 suc h that 9 < — and 

supprj = {re 10 : n ^ r ^ r 2 , -9 < 9 ^ 9 }; (4.21) 

(hi) for every nonnegative integer J, we have a nonstationary tight wavelet frame in /^(R 2 ).' 

E K/-WI 2 + EEE \(f^ij\ 2 = WfWUu*), V/ G L 2 (M 2 ), (4.22) 
kgz 2 j=J l=\ kez 2 

where Sj := m2™ ; ['J ^ e floor function, and all are real-valued functions in the Schwarz 
class satisfying the following properties: 

(1) ipi' is defined as follows: for r ^ and 9 G [— 7T, 71"), 

- j^e^ + rji-re^), if 9 G [-2"^, 2-^tt), 
^ } \0, i/0 G [-7r, 7 r)\[-2-L«J vr ,2-L«j 7r); ^.23) 

(2) all other generators ijp' are obtained via rotations from ip^°: for all £ = 1, . . . , Sj, 

^ e (re i0 ) := ^'°(reV L/ " J,r( '- 1)/ro ), r ^ 0, 9 G [-vr, tt); (4.24) 
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(3) the support of ip^f^-k ^ as ^ wo P ar ^ s which are symmetric about the origin and each part 
obeys width ~ length 1 ^ p . More precisely, for all k G Z 2 and j ^ ; 

supp^^ k = {re ie , -re ie : 7Pr x ^ r ^ 2 J 'r 2 , -2^ P ^6 Q <: 6 ^ 2 j '-LpjJ 6> } ; (4.25) 

(4) are compactly supported functions in C°°(]R 2 ) vanishing in a neighborhood of the origin 
and there exist 27rZ 2 -periodic trigonometric functions hj^ in C°°(T 2 ) such that 

5^(20 = KiiOkO, (GR^ = l,..., SjlJ eNU {0}. (4.26) 

Proof. To prove Theorem [TU we construct some special 27r-periodic functions. For a positive 
integer n and any < e ^ j-, by |IE1 Lemma 4.2 or (4.2)], we can construct a 27r-periodic 
nonnegative function a n ^ £ in C°°(T) such that 

«n, £ (0 = l, Ceh^+^^-s] and 0^(0 = 0, £g [_7r,7r)\[-^-e,^+e] (4.27) 
and 

2n-l 

EMH!)l 2 = i, eet. (4.28) 

Let (p and 1/? be constructed in the proof of Theorem [2] with d = 1 and M = 2. Or equivalently, 
by [HI Lemma 4.2 or (4.2)], we can easily construct two real- valued even functions ip and i/> in 
C7°°(R) such that their supports are contained inside [— tt, tt], tp takes value one in a neighborhood 
of 0, ( 14.14p holds with N = 1/2 and d = 1, and there exist 27r-periodic trigonometric functions a 
and b in C°°(T) such that 

0(2£) = a(O0(£) and ^(2£) = *>(£)<£(£), ^M. (4.29) 

Now we define 

<p(re ie ) = <p(r), i>(re ie ) = ^(r) and 77^) := ^(r)a m , e (0), r ^ 0, G [-tt, tt). (4.30) 



One can easily check that the supports of all ipi' 1 are contained inside {—tt, Tr) d and ( 14.201) holds 
with d = 2 and N = |/ 2 for j G N U {0}. Define 27rZ 2 -periodic functions a, bj j0 on [—tt, 7r) 2 by 




For £ = 1, ... , S j-, we define b j>£ (re ie ) := b jt ( re i ( fl - 2Lw ' Jfl ' r (*- 1 )/" 1 )) for < r < tt and b ij£ (re ie ) := for 
otherwise re* 6 G [— 7r,7r) 2 but r ^ tt. By Corollary [17] using a similar proof as that of Theorem [2j 
we see that all the claims can be easily checked. □ 
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